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ELECTROSTATIC  ION  INSTABILITIES  IN  THE  PRESENCE  OF 
PARALLEL  CURRENTS  AND  TRANSVERSE  ELECTRIC  FIELDS 
1 .  INTRODUCTION 

It  is  veil  known  that  an  equilibrium  current  parallel  to  the  magnetic 
field  may  lead  to  the  excitation  of  obliquely  propagating  electrostatic 
ion  modes  in  a  plasma,  such  as  the  ion-cyclotron  (EIC)  mode  in  the 
collisionless  (Drummond  and  Rosenbluth,  1962;  Kindel  and  Kennel,  1971)  and 
in  collisional  limits  (Milic,  1972;  Chaturvedi  and  Kav,  1975)  and  ion 
acoustic  waves  (Chaturvedi  et  al.,  1987).  These  instabilities  have  been 
extensively  studied  in  the  context  of  the  auroral  ionosphere  (Kindel  and 
Kennel,  1971;  D'Angelo,  1973;  Chaturvedi,  3976;  Sa tyanarayana  et  al . , 


19E5;  Providakes  et  al . , 


Cha’.u:  ved;  e  t  al.,  198' 


Radar 


observations  (Fejer  et  al.,  1984;  Haldoupis  et  al.,  1985;  Villain  et  al . , 
1987)  and  in-situ  rocket  measurements  (Kelley  et  al.,  1975;  Ogava  et  al . , 
1981;  Yau  et  al.,  1983;  Bering,  1984)  have  detected  these  waves  in  the 
auroral  ionosphere.  It  has  been  shown  by  Kindel  and  Kennel  (1971)  that 
the  EIC  instability  has  the  lowest  threshold  of  excitation  among  the 
various  current  driven  ion  instabilities  for  the  topside  ionosphere. 
Thus,  these  observations  have  been  interpreted  as  the  modes  excited  by 
field-aligned  currents  (FAC).  Often  however,  along  with  the  FAC, 
transverse  d.c.  electric  field?  have  also  been  observed  (Fejer  et  al., 
1984;  Providakes  et  al.,  1985J,  and  the  threshold  values  for  excitation  by 
FAC  alone  are  likely  to  be  exceeded  only  during  periods  of  strong 
geomagnetic  activity  (Providakes  et  al . .  1985;  Chaturvedi  et  al . ,  198?  )  . 

Recently,  it  has  been  demonstrate  that  the  presence  of  transverse 


localised  electric  (TLE)  fields 


result  in  the  excitation  of 


electrostatic  ion  modes  (Ganguli  et  al-.  l°85a.  1985b. 1988).  These  mode? 

have  recently  been  seen  in  PIC  simulations  conducted  by  Nishikava  e'  al . 

[198"(.  This  nonlocal  instability  is  a  result  of  coupling  between  the 

negative-energy  wave  in  the  localized  electric  field  region  with  the 

positive  energy  wave  outside  it.  These  studies  were  for  collisionless 

plasmas  and  largely  for  k  >>  k  waves:  but  the  important  role  of  the  TLE 

fields  in  a  realistic  study  of  the  low- f requency  ion-wave  excitation  in 

the  auroral  ionosphere  was  demonstrated.  In  this  paper,  we  study  the 

combined  effects  of  a  TLE  field  and  an  FAC  on  the  stability  of  the  ion- 

modes  .  Ue  also  include  the  electron  neutral  ( v  )  and  the  ion  neutral 

en 

(v.  )  collisions.  Ve  find  that  the  inclusion  of  TLE  field  can  have 
in 

notable  effects  on  the  dispersive  characteristics  of  the  ion-modes.  For 
example,  the  real  frequency  of  the  unstable  modes  can  have  values, 
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depending  on  parameters  such  as  the  TLE  field  strength,  that  vary  from 
below  Sj  to  values  greater  than  0^ ,  the  ion  cyclotron  frequency.  Hore 
importantly,  we  find  that  the  mode  excitation  .flay  occur  for  sub-threshold 
parallel  current  (electron  drift)  values  when  the  TLE  field  is  included. 
These  results  may  be  of  interest  for  the  ion  wave  observations  in  auroral 
ionosphere,  especially  for  those  observations  where  it  appears  that  the 
ambient  conditions  set  the  threshold  value  of  the  parallel  current  for 
excitation  of  the  ion  waves  above  the  level  that  one  might  expect  to  be 
attainable  under  normal  c i r cums tances .  A  combination  of  a  parallel 
cut  rent  ar.d  a  TLF.  field  car.  lead  to  exc  i  r  a :  :  on?  of  the  ton  waver  unde: 
much  less  severe  conditions. 

2.  THEORY 

Ve  wish  to  study  the  stability  properties  of  a  weakly  ionized 
collisional  plasma  which  is  typical  of  the  auroral  ionosphere, 
characterized  by  an  equilibrium  current  flowing  parallel  to  the  ambient 
magnetic  field  and  a  localized  transverse  d.c.  electric  field  (see  figure 
1).  Ve  describe  the  particle  dynamics  by  using  the  Vlasov  equation  with  a 
BGK-type  of  model  relaxation  collision  term.  The  equilibrium  current  is 
predominantly  carried  by  thermal  electrons  drifting  along  the  magnetic 
field.  BqZ,  with  a  drift  velocity,  V^z.  The  transverse  d.c.  electric 
field,  Epfx),  is  localized  along  the  x-axis  to  a  distance  and  is 
piecewise  continuous.  Thus,  the  ions  and  electrons  both  have  an 
equilibrium  cross-field  drift  (V^y,  V(„  =  -cE^  B-).  which  is  also  localized 
in  the  x-direction  to  the  distance  L_  .  In  reality,  however,  the  profile 
of  Vq  is  no  longer  of  the  square  hat  shape  and  is  rather  smooth  due  to  the 
f  ini  te-Lamor-radius  (FLR)  effects  at  the  boundaries  of  the  electric  field 
region.  For  simplicity  we  shall  use  a  sharp  profile  for  in  the  present 
study,  and  shall  subsequently  relax  this  condition  in  the  future.  A 
similar  approach  was  taken  in  developing  the  col  1 isionless  TLE  only  theory 
(Ganguli  et  al.  1985a,  1985b,  1988],  where  it  was  found  that  the 
transition  from  the  ideal i zed  square  hat  profile  to  a  more  self  consistent 
smooth  profile  resulted  in  quantitative  but  not  qualitative  changes  in  the 
theoretical  predictions.  The  effects  associated  with  the  Hall  and  the 
Pederson  mobilities  are  small  for  weak  collisions,  especially  when  the 
velocity  shear  is  weak;  and  are  therefore  neglected.  Our  treatment  is 
general  in  as  much  as  it  includes  collisions  and  thus  is  appropria'e  for 


2 


lov-altitude  applications;  and,  at  the  same  time,  the  collisionless 
results  (applicable  at  higher  altitudes)  can  be  readily  obtained  by 
setting  the  collision  frequencies  to  zero. 

The  dispersion  properties  of  the  collisional  ion  cyclotron  modes  in 
the  local  limit  and  in  the  absence  of  a  d.c.  electric  field  are  described 
by  the  relation 

1  +  e.  4  t  =  0  (1) 

i  e 

where  expression  for  e .  ( j  =  e.i)  is  given  as  jClemmov  arid  Dougherty.  1  *>69 ; 

J 

Kindel  and  Kennel,  1971) 


e .  = 


*  *  (VkllV.j)Z[(“j  ♦  ngj) 


,  |/k..vt. 


OD 

1  +  ifv./k  v  1  I  T  (u.)Z  (go.  +  nC.]/k..v..l 

V  1  h  tjJ  „  nv  j ;  LI  J  jJ  II  t  j  J 


n  =  -® 


(2) 


The  notations  used  here  are  standard: 
number;  go 


n  refers  to  the  harmonic 
go..  *  iy; 


go  *  iv  where  go  is  the  complex  eigenf  requency ,  go 

J  2  1*2  r 

Vi  =  <T  ./4II  n j e  )  is  the  Debye  length  of  the  species;  and  and 

m.  are  respectively  temperature  (in  energy  units),  number  density  and  mass 
J  1/2 


of  the  species  j ; 


J 


(2T./m. ) 
J  J 


is  the  species  thermal  velocity;  \k 


their  collision  frequency  with  neutrals:  p.  =  v./2.,  the  gvroradius;  T  = 

J  J  J  tr 

I  exp  ( -b .  ) .  where  I  is  the  modified  Bessel  function  of  order  n;  b.  = 

n:  j  J  n  J 

/ 2:  and  2(C)  is  the  plasma  dispersion  function  f3ook,  1986]. 

Ve  now  introduce  a  localised  transverse  d.c.  electric  field  (see 

figure  1).  This  modifies  the  dispersion  relation.  Region  I  (x  <L^/2), 


over  whicn  the  electric  field  :s  localised. 


V=voy:  “i =  * 1  v 

=0;  go.  =  go  *  i  v.  ;  go 
Oi  ’  i  l  e 

Ve  shall  assume  thai 


u)  -  GO.  - 

e  1 
=  go  +  i  u 


k  7  , . 
z  d 


is  characterized  by  Eq  =  Eq.x; 
In  region  II  we  have  En=Q; 


k  7,.  Here  go.,  =  go 
z  d  1 


Vo- 


the  density  gradient  is  negligible.  Strictly 
speaking,  the  parallel  current  is  also  of  finite  width  in  the  x-direction. 
The  typical  scalelength  of  this  width,  .  is  assumed  to  be  either  Lc  >  Ly 
or  -  Ly .  Ve  discuss  implications  of  these  inequalities  in  the 
discussion  section.  For  a  more  rigorous  description  of  the  equilibrium  we 
refer  to  Ganguli  et  al.  (1988). 


Ve  follov  the  approach  adopted  by  Ganguli  et  al.  [ 1 985a )  in  deriving 
the  nonlocal  dispersion  relation.  Ve  use  the  quasi-neutrality  assumption 
(c^  *  t.  ;  0).  As  a  result  of  the  x-dependence  of  the  equilibrium  the 
perturbed  quantities  cannot  be  Fourier  transformed  in  x.  Instead  ve 
obtain  a  differential  equation  in  x.  For  the  specific  profile  of  the 
electric  field  under  consideration  (see  figure  1),  the  modes  are  described 

by 


,  for  Region  I 


where  £,  = 


-  Kj  j  *j(l)  =  0 


fe  *  hi)  '  0 

o  c. 


x/p.  and  kt  =  QT/AT  and 
l  I  II 


,  for  Region  II 


„I  •  c  1  .  t  rn  (; 

V  n 


ja.  *  nC. 

/k  v  .  Z  [-7 - - 

"  1  knvi 


T  D  1  -  i 


,  .  a>.  *  n2. 

v.  /k..v.)  r  r  z  p-r - - 

i  II  l ;  n  v  k,,v. 


.  f  f  -  .  A'.  *  nB.  .  1 

1  U  -  | .  ,  ■  l  l  > 

L  =  -  ^  jCi.w  .  IV.  k,,v  - : - |> 

1  c  I  n  1  1  tc  v  .  ;  | 

L  ji  i  i 


nl,  r,  [■  ,,  f^i  ~  °^i j 

*  T  D':  Fr,  !1Vi  k,Vi)Z'  k  v.  I 


C  =  | 1  *  i  fv  /k*v  )z (w  zk  v 
I  v  e  f  i  e  /  \  e  e 


D  =  <1 


u>  /k..v  |Z  fu>  k  -j  ) 

e  ne,'  \  e  1 1  e , 


Only  the  n  =  0  harmonic  for  the  electrons  is  considered  and  T  '  =  8T  /8b. 
1  n  n 


;  . . . ' 


In  the  above,  T  =  T./Te,  and  Kjj  is  obtained  by  setting  Vp  «  0  in  tj. 
The  nonlocal  dispersion  relation  for  the  even  modes  obtained  by  matching 
the  logarithmic  derivatives  of  the  solutions  of  (4)  and  (5)  at  x«Lv/2  (for 
details  see  Ganguli  et  al.  1 1985a | )  is  given  by, 

-  Kj  tan  (Kj/2e)  =  i  tCjj  (9) 

where  e  =  p./L  .  A  similar  relation  can  also  be  obtained  for  the  odd 
1  v 

modes . 

The  well-known  limits  of  (9)  are  straightforward  to  obtain.  In  the 
local  approximation,  and  for  =  0,  we  recover  the  current  driven  EIC 

instability  in  the  collisionless  (\m=0)  and  the  collisional  (v^*0)  limits. 
Further,  for  =  0  and  =  0,  (9)  reduces  to  the  dispersion  relation 

given  by  Ganguli  et  al.  (1985a]  to  yield  an  electrostatic  ion-instability. 

Here  we  study  the  combined  effects  of  the  two  aforementioned 
processes,  i.e.,the  simultaneous  presence  of  the  FAC  and  the  TLE  field. 
There  is  evidence  that  except  for  periods  of  moderate  geomagnetic 
activity,  the  observed  values  of  the  field-aligned  currents  fall  in  the 
sub- threshold  domain  for  the  collisional  EIC  and  IA  instabilities  (Fejer 
et  al.,  1984;  Chaturvedi  et  al.,  1987).  Also,  there  are  several 
suggestions  that  the  observations  indicate  the  presence  of  transverse 
electric  fields  in  the  auroral  regions  [Fejer  et  al.,  1984:  Basu  et  al., 
1984;  1985;  Prikryl  et  al.,  1986].  Ve  therefore  investigate  the  role  of 
the  TLE  field  on  the  current  driven  ion  modes. 

3.  NUMERICAL  RESULTS 

We  proceed  to  evaluate  the  nonlocal  dispersion  relation  (9).  We 

consider  ion  and  elecrron  tempera -.ures  to  be  equal  (t  =  i).  Other  typical 

parameters  used  are  v  =  v  /2  =  12.  v  -  v./2.  =  0.033,  £  =  0.1  and  u  = 

K  e  e  i  ill 

k  ,/k  =  0.15.  We  studv  an  0*  ior.-plasma  svstem  i.e.,  u  =  m./m  =  29392. 
i|  *  r  l  e 

Figure  (2)  shovs  the  real  and  the  imaginary  parts  of  the  frequency 

O  __  _ 

(u  and  y)  as  a  function  of  b(  =  k  ^  p . “  Z )  for  VA=V„/v.=0  and  V  ,  =  V,/v.  =  30. 
r  ±  l  0  0  l  d  d  l 

The  mode  is  stable  For  the  k(  domain  of  interest.  The  real  frequency 
shows  the  expected  weak -dependence  on  b.  This  is  in  agreement  with  the 
previous  studies  which  suggest  a  higher  threshold  value  of  for  EIC  wave 
excitafion  ( Sa tyanarayana  et  al.  [l??c'],  Providakes  et  al  [1985]). 


Next  we  consider  V,  *  0  and  let  9-..  -2.8,  v  «  0,  v,  ■  0.  From 

a  U  e  1 

figure  (3)  ve  see  that  the  real  frequency  of  the  mode  has  a  roughly  linear 
dependence  on  b.  Ve  find  that  the  mode  is  stable  for  the  b-domain  which 
we  have  scanned  and  which  is  of  interest  to  us.  Thus,  for  9,  =  30,  9A  =  - 

_  _  <3  U 

2.8,  u  =  0.15,  =  12,  v.  =  0.033,  the  nonlocal  electrostatic  instability 

discussed  by  Ganguli  et  al.  (1985)  and  the  collisional  current  driven  E1C 

instability  are  both  stable  individually. 

Nov  ve  combine  the  tvo  drifts  (9^  and  9^)  for  the  values  given  above 

and  solve  (9)  numerically.  The  results  ( w  and  y  plotted  vs.  b)  are  shovn 

in  figure  4.  A  comparison  vith  figs.  (2)  -  (3)  reveals  that  the  mode 

frequency  follovs  the  approximate  pattern  of  the  nonlocal  electrostatic 

instability  (fig.  3)  but  the  grovth  rate  is  positive.  Thus,  ve  find  that 

in  presence  of  a  parallel  electron  drift  ( ) ,  the  nonlocal  electrostatic 

instability  exhibits  grovth,  even  in  the  presence  of  collisions.  Based 

upon  this  result,  ve  suggest  that  in  the  collisional  lov-altitude  auroral 

ionosphere,  the  observed  lov  frequency  plasma  vaves  may  have  been  excited 

due  to  the  simultaneous  presence  of  both  the  transverse  electric  field  and 

the  parallel  current.  A  similar  conclusion  for  the  collisionless  (high- 

altitude)  case  can  also  be  dravn,  indicating  that  excitation  of  the  ion- 

modes  jointly  by  the  parallel  current  and  the  transverse  electric  field  is 

possible  for  sub- threshold  electron  current. 

Ve  find  that  the  effects  of  finite  channel -vid th  (L  )  on  the  nonlocal 

c 

ion  instability  are  relatively  small  if  >  L  .  In  figs.  (3)  -  (4),  ve 
have  assumed  >  Ly.  For  corner  :s^r .  ve  plot  y  vs.  b  in  fig-  u  for  the 

case  =  f...  (dashed  curve).  Ve  see  that  the  grovth  rate  in  this  case  is 
of  the  same  order  as  before  (I  >  L,  )  though  slightly  larger. 
Qualitatively  the  reason  for  the  relatively  smaller  influence  of  finite  L 

c 

on  the  nonlocal  ion  instability  in  cases  vhen  L  >  L  is  that  the  va ve 

c  -  v 

packet  (localized  around  the  electric  field)  samples  the  region  of  near 
peak  current  in  space. 

The  value  of  the  |Vq|  =  2.8  used  in  the  numerical  computations  above 

is  larger  than  the  typical  observed  value  of  j9^j  -  0.5.  In  figure  (5)  ve 

nov  consider  9»  =  -0.5  and  V  =  30  and  25  vith  u  =  0.15  and  0.17. 
0  d 

respectively.  The  rest  cf  the  parameters  remain  unchanged  from  the  figure 
(4).  Nov  ve  see  a  rather  coherent  electrostatic  vave  grovth  around  the 
ion  cyclotron  frequency  for  sub-critical  values  of  9^. 


4.  DISCUSSION 


Ve  have  shovn  in  this  article  that  the  simultaneous  existence  of  a 
localized  transverse  electric  field  and  a  parallel  current  can  excite  the 
oblique  ion  modes  even  though  the  electric  field  and  the  current  values 
are  separately  stable.  The  parameters  ve  considered  are  typical  to  the 
lov-altitude  auroral  ionosphere  (the  upper  E-region).  In  this  region 
various  experiments  have  indicated  observations  of  EIC  or  EIC-like  ion- 
modes.  These  observations  are  frequently  attributed  to  a  field-aligned 
cur ren t -dr i ven  EIC  generation,  but  often  the  current  is  below  the 
threshold  value  requited  for  excitation  |ree.  e.g.  ,  Fejei  et  al.,  198m], 
And  the  frequency  of  the  modes  may  also  display  variations;  it  may  be 
equal  to  the  ion-cyclotron  frequency  or  sometimes  different  from  it 
(Prikryl  et  al.,  1987 1 .  In  absence  of  detailed  data  on  the  ambient 
parameters,  it  is  not  possible  to  obtain  close  comparisons  with  the 
observations.  However,  for  the  parameters  we  have  considered,  it  appears 
that  the  excitation  of  the  transverse  electrostatic  modes  by  a  combination 

of  a  TLE  field  and  a  parallel  FAC  can  be  a  possibility.  The  parallel 

2 

currents  (corresponding  to  tens  of  uA/m  ,  i.e.,  -  30)  and  the 
transverse  electric  fields  (corresponding  to  100  mv/m,  i.e.,  |Vq|  -  0.5) 
correspond  to  the  typical  values  that  are  believed  to  have  been  observed 
in  the  regions  of  auroral  wave  activity.  Further,  the  range  of  unstable 
mode  frequencies  varies  from  below  2.  to  values  greater  than  2.  and  may 
have  relevance  to  the  auroral  situations  where  the  observed  wave 
frequencies  sometimes  are  different  from  the  ion-cyclotron  frequency. 

The  mechanism  of  mode  excitation  in  the  combined  presence  of  and 
V  is  of  interes*  to  other  situations  as  veil.  For  example,  a'  higher 
auroral  altitudes,  the  observations  of  EIC-like  modes  (and  lover  hybrid 
(LH)-like  modes)  have  been  made  where  the  TLE  fields  may  have  been  non- 
negligible  (Kintner  et  al.,  1979].  Ve  are  presently  conducting  a 
parameter  study  for  this  case.  In  addition,  the  profile  of  electric  field 
chosen  by  us  in  this  paper  is  an  idealization.  Preliminary  i nves t i ga t i ons 
with  a  smooth  profile  suggest  that  our  results  are  not  severely  modified. 
Detailed  results  vill  be  provided  in  a  forthcoming  publication. 

In  conclusion,  we  have  shown  that  the  sub- threshold  parallel  currents 
and  stable  transverse  electric  fields  may  combine  to  result  in  the  growth 
of  elec t ros ta t i c  ion  modes  in  a  collisional  plasma.  Ve  have  applied  our 
results  to  the  lower  altitude  auroral  ionosphere  and  find  that  reasonable 


values  of  parallel  currents  and  transverse  fields  may  result  in  the  ion 
■ode  excitation  under  relatively  less  severe  conditions  of  geomagnetic 
activity.  Detailed  results  including  the  collisionless  case  applicable  to 
higher  altitudes  will  be  presented  in  a  future  publication. 
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